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ABSTRACT ■ 


The problem of potential steady subsonic flow around 
complex configurations is considered. This problem requires 
the solution of an integral equation relating the values of 
the potential on The surface of the body to the values of the 
normal derivative, which is known from the boundary conditions. 
The surface of the body is divided into small (hyperboloidal 
quadrilateral) surface elements, 2^ , which are described in 
terms of the Cartesian components of the four corner points. 

The values of the potential (and its normal derivative) within 
each element is assumed to be constant and equal to its value 
at the centroid of the element. This yields a set of linear 
algebraic equations. The coefficients of the equation are given 
by source and doublet integrals over the surface elements, X ; . 
Closed form evaluations of the integrals are presented. 
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SECTION I 

FORMULATION OF THE PROBLEM 

1.1 Introduction 

A general theory for compressible unsteady potential aero- 
dynamic flow around lifting bodies having arbitrary shapes 
and motions is given in Refs. 1 and 2. Application to finite- 
thickness steady and oscillating wings in subsonic flow is 
given in Refs. 3, 4 and 5. Here, a general numerical formu- 
lation for complex configuration in steady subsonic flow is 

considered. By using the Prandtl-Glauert transformation, the 

1 2 

incompressible flow is obtained. ' Hence, for simplicity, 
only the incompressible flow is considered here. In this 
case, the problem is governed by the Laplace equation with 
prescribed normal derivative on the body (exterior Newman 
problem for the Laplace equation) with an additional compli- 
cation due to the presence of the wake (of unknown geometry). 

The method is described with the emphasis on the aerodynamic 
applications, but it is applicable to different physical pro- 
blems as well (see Subsection 6.3). 

The problem of the evaluation of the steady, incompressible 
potential aerodynamic flow around an aircraft of arbitrary 
configuration can be analyzed by solving the integral equation 

(1 - 1) 

2 5 

where f is a surface surrounding the aircraft and the wake. ' 
Subsonic oscillatory flow is considered in Appendix C. 
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For the moment, it will be assumed that the wake does not 

exist. The effect of the wake is considered in Section 5. 

The value of ~~~ is obtained from the boundary condition 
o n 

(tangency condition) 


3 $ c, 

On “ ( f) ' Uco ( H - 


) -o 


( 1 . 2 ) 


29 _ _ 

an ~ ^ ~ H ♦ t 


(1.3) 


The integral equation can be studied by dividing the 
surface <5* into N small finite elements <5~; to yield 


l, ( ) f” -1 T * 


(1.4) 


Applying the mean value theorem one obtains 




'* *Vi~ Sl % JJ n ’ y 


(1.5) 


where ^ is a suitable mean value of <j> inside the element , 

which will be approximated by the value of (j> at the centroid 
P v of the element, <$ ^ 

By satisfying Eq. (1.5) at the centroid, Q f the 

element <51 , (h = 1, 2....N) yields 


*-^4, IJ 


+ L % f/”’ b, j dS * ( h ‘ t 2. i. M) 


(1.6) 
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where r^ is the distance of the centroid of the element 
from the dummy point of integration in the element <? 
Equation (1.6) is equivalent to* 


where 


and 


with 


[ W- = \K) 

<f* * 


(1.7) 


( 1 . 8 ) 


(1.9) 



1.2 Surface Geometry 

Let the geometry of the element 6* K be described by 

p > P -5') <1 ’ 11 

where 1 and v are the generalized curvilinear coordinate. 

Then the two base vectors ft. are given by (Fig. 1) 

— ■> 

- P 

a. = ■ — — . d.i2 

The unit normal to the surface is given by 

ft. * G* (1.13 

a, 11 0. 


The effect of the wake is not considered here (see Sections 


5 and 6) . 
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and is directed according to the right-hand rule (Fig. 1) . 
The surface element is given by (Fig. 1) 


d<$ 


(X, d-J 1 x 


| d-f ’d-f 2 


(1.14) 


1.3 Expressions for b^ and c 


'hk 


Combining Eqs . (1.10), (1.13) and (1.14) yields 


bh* = j- || (a, x I) — d-f 

Similarly, combining Eqs. (1.8), (1.13) and (1.14) yields 


(1.15) 


where 


= 2TT ([ a ' * ‘ ^ (y ) d-f 'cj-f " 




-40 


Q, * Qx ■ r A 


df'd-f 




( 1 . 16 ) 


(1.17) 


* -X, 

J - »K 
2 - 

In Section 2, these expressions are evaluated under the 
hypothesis that the surface element is a portion of a hyper- 


boloid. 
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SECTION II 

HYPERBOLOIDAL ELEMENT 


2.1 Introduction 

Consider the equations 

"X - ■'X c + 'X,-f + y %i r ) + ^ ^ 9 

^ "f » 'L 1 . ^ -f 7 

J ’ K * *■-? - M ♦ j.-f, 

or, in vector notations 


(2.1) 


p ■ P. ’ M * p>*7 - p,-?7 (2 - 2) 

This represents a hyperboloid. The lines r J = const and = 
const are clearly straight lines. Consider the hyperboloidal 
element defined by the above equation with 

■I <-1 < I 
- 1 17 - I 

(2.3) 

The centroid of the element is p c ( 0). The corner 

points of this element are 

P« * P< * P, * p.- f, (-?■-!, 7--.,) 

P.- = Pc - P, - P. - P> (-f = . 1, 7 = - 1) 

P-. ' Pc - P, •* P.-p, (-? = -! , 7-') 

P- * Pc - P, - P. - P, 


(2.4) 
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The inverse relation is 

?<. - i (P..- ?..* j>,- P..) 

?■ ■ + i- - - r.) 

?, . i( i-i - f„- i) 

(2 • o ) 

Note that the four boundaries of the element ( - t / } + | ) 

are straight lines given by 

?* <P.*f>)7 -I 4 ? * | 

P-- ( pc - P.) * <?>-?,)? -I 4 •? £ I 

P= <?.* + -I si i i 

p-- (pc- p,)^ -i sit i 

( 2 . 6 ) 

Next, assume that the surface of the aircraft is divided 
into curved quadrilateral elements with four corner points 
p ++ , P+_' P_ + > p__. Then, as mentioned in Section I, these 
elements can be replaced by the hyperboloidal element (des- 
cribed above) which goes through the four corner points p ++ , 

p + _ , P_ + / P (see Fig. 2). It may be noted that the surface 

is continuous since adjacent elements have in common the 
straight line connecting the two common corner points. It 



7 


may be noted also the p^, is the centroid of the hyperboloidal 
element 6"^ and hence it will be indicated as 

p , p (2.6) 


2.2 Geometry of Hyperboloid Element 

The geometric quantities introduced in Section I can be 
written for the hyperboloid element described above. Letting 

i' * 1 , < * ? (2.7) 

Equation (1.12) yields 

a, = — = p, ? 3 <J (2.8-a) 

-* 7, tf -* _ 

= jf - ?, + (2.8-b) 

This yields 

o. * «. = ( p, « ?>?)< <p. « p,-j) 

“ P> ‘ f> * Pi * P> - p, « R n 

' (2.9) 

since x = 0. In components notations 



I 1 

j 

£ 


\ l 

j 

i ' 

J 

K 

j 

* 

n 0, - “Pet 


f. 

*ii 

+ Det 1 

' X, 


1 1 

Pet | 

X J 



« 

( 

. ** 


iJ 

l 

*9 

* 

hi 

? , 




| W.5.- 'j 


( 2 . 10 ) 
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In particular, for the first component, one obtains 


= <2 ' L1> 

Note that, with present notations 

<L - f £ -«*) f, ( - ih) - 

?--h*p-P »p -p -p,i - P^-p^f 

" P- * Pt - P,9 i- p, (2.12) 


where 

P. . p‘*’_ p <M 


(2.13) 


is the vector connecting the centroid p ^ of the element <S h , 
to the one, p , of the element <5^ . Hence 


^ * fl, X 

(p, » p. i- p, * p»-f * %> f> i )- ( p. <• p.-f.p. ?j-? ? ) 


= ( p. • p, < p.) * <p. • P. * %)i - ( p.- ?, * pj 7 

* ^ > ?.>? ? + (p»-p, >p,)t ?-(p> •?.«?.>?? 


since 


(2.14) 


?r Pi x P» j ?i • P' * P» * P- ' ?, « ?, -- ?,• f, * f, 

- viw, * Ps ■ P> " P. - O (2 . 15 > 

Noting that 

= f* * ?, * r - ?,* h - - £ ? 3 


(2.16) 
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Equation (2.14) reduces to 

}■ S,*o. . fp.-p,.p.)* (P.- P,*bH 
+ (p. • p, « p»);- ( p, • ?.« 


In component notations, Eq . (2.17) reduces to 

r ~ - ^ £» 7( o 

^ • a,« - Dtt ^ ^ ^ + Det „ 

■x . (.(. 


T, i. 

+ Pet , ^ 

•X, *. }, > 

"X i ^ > 5j| 


^ ^ lo [ 

*. 1. Si 

^ h n - Vet 


0v 7 

1 j 


"5 7-. 

" I 


% i. 
h h 
p. 

i> h 


Finally, according to Eq. (2.12) 


^ = I n | * i j | -|p.- p.-j -r p, 7 , 

-{ P. * P. ’ P, • Pi P.‘ P, j‘t P 4 • P ( -J ‘7* 

+ iCp.-p,t * P. ■ p. 7-. p, • p.)f f,-P,l‘7*F,-f,1p]} 


or, in components notations 
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2.3 Expression for and 


By combining Eq. (1.15) with Eqs. (2.11) and (2.20), one 
obtains 


i t 


b 


hh - + 


2-TT 


( lr, ty,)* ($,£} ~ lr, \ 


•1 'I 


* t + ($,* * y>) r ty-f •>) + (} 9 t t fof?) J x <J-f Ay 

( 2 . 21 ) 

Similarly, combining Eq. (1.16) with Eqs. (2.18) and (2.20) 
yields , 


Ch * = ‘ rR j) ^ et 

-i -i 


X, 

1* 


Xi 


£. 

X. 



X J 

h 

S-J 

*»-??/♦ (S'. 

- 


(2.22) 


The integration of Eqs. (2.21) and (2.22) is discussed in 
Section 3 and 4, respectively. 
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SECTION 3 


DOUBLET INTEGRAL 


3,1 Integration with respect to > 

Consider Eq. (2.22), which can be rewritten as 


hi k — — 


fV . 

) ) ir » ^ 


De-t 


*• h 

Xl K 

* > $-> 


- De-t 


^ 


(3.2) 


^ I = I 3. + f i, | = (&..+ 2 ( 3 '3) 


with 


= 2-i • £; 


(3.4) 


where 


P. * ? ?, 

P, * ? P, 


(3.5) 



12 


The indefinite integral with respect to in Eq. (3.1) is 

given by 





(3.6) 


In order to obtain the relation between n^, and 
consider the derivative of the expression in bracket in Eq. 
(3.6) and equate it to the integrand of Eq . (3.1). This 


yields 


3 

55 


n»+ n,5 

. r h 




= [ h| (r h ( Zot + Zn "f )] 

- ~ 3 [ n. ( a,. + %,-% ) - n. C2., "1)J 

• h 



-f ( W.- m.f) 

'h 


This implies 



1 

° c ^t 

< 

C 

*»»— — 

•» 

>n. ) 


2,, , 


n 0 , 


hi , , 


or 





£ 


0 I 









(3.7) 


(3.8) 

(3.9) 
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Note that (see Eq. B.l) 


V2, 2.- 2, 

£,• i. a, -a, 


D - Detj 2 ” ’ ,2 "' = D„[ V2 * 2 ’"^ 

2„J li,' i. a.-i., 

- I 2.« i, | j 

Finally, Eq. (3.9) yields 

110 = + 

■ ^ +■ "f a,)* a.w,) 

; p H 

since, according to Eqs. (2.12) and (3.5) 

} - P . + -?P, *U, --jyp, 

- 2, -• ■? 2, ' 

and, having defined 

Q = a, 1- a, m , 


2. 7 

X, 

h 

2,7 

A> 

i 

- 2. 

X, 


-2, 

Xj 

h 



(3.10) 


(3.11) 


(3.12) 


0 W** 

0_ X l + 7X } 

-2. K 

~2| X} ^3 ti 
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1 X „ + X >- 



- JL 

X. 

i* 

it 


l 

h 



x„ 


*. ' 

+ 3, j 

X, r ^ Xv 




s X i 


K J 


= - a. (a.- p, » p,)+ 2, (f.- a,* pj 

s - 2. ( 2. • P,» p,) - 2,(2, • R « pj 

= - ( p.^7f.)[ ( P."?P.)' P, « pj 
- ( p, - 1 ?,)( (FrlF,)- P.J 


(3.13) 

Finally, combining Eqs . (3.6), (3.11) and (3.13) yields 

l v (n) = _ p.«pj J 

- 2, 2 IS.+fs.) 


(3.14) 
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an alternative form is (note that * P, * P i = S 0 * x P 3 
and £, • R x R » 2 ,* * P, ) 

• , , = [( 3,1- fa.) p, - (a,,.faj pj ( a. « i.) 

" ~ 2 ... 2„- ' | 2.1- fi.l 

= (a. "fa,) - (a, p, - a, p,)-(i. >i,) 

|2."5l,| (a.*£,K£.«2,) 


I |l(a„< 5,| : 


(f a,) 

■s 


If I if 




i 

Iff If Afl,|* 


[a* ( f.* f)(f * 

(3.15) 


3 ♦ 2 Procedure to be Avoided 

It may be noted that, according to Eqs. (3.1) and (3.14), 

C, is given by the sum of two integrals with respect to rj is 
hk 
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of the type 


1 _ ( ! d n 

15 “ J f / 


(3.16) 


where 

while 


M}(y) = I is a polynomial of third degree in ^ 


N«<?> * $ 7* - 1 a. - a. | >o 


(3.17) 


is a nonnegative polynomial of fourth degree. If N 4 ( ) were 
always positive, then the rational function could be replaced 
by a polinomial 


M» {<?) 1 

— ^ . c. + (3.i8) 

KM# 7 7 

and the integral in Eq. (3.16) could be easily evaluated by 
using the recurrent formula (Eq. A. 6) 




+^7^ r 

(3.19) 


together with (Eqs. A. 9 and A. 10) 






(3.20) 
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However, if for a certain value, y* , of *) the denominator 

is equal to zero, Eq. (3.18) cannot be used. It may be noted 

that this implies at ? * 9+ f 1 1, *2L,\*0 or, that 2 , * P* + p, 

is parallel to ^ f 4 . As is evident for Fig. 2, this is 

Ph (*) 

the case if the poinybelongs to the line y s rj * . It may 

be noted that, as shown by Eq . (3.15), the numerator has a 

single root at J - y while the denominator has a double root 
at y 1 . In order to avoid this problem, it is convenient 
to follow the procedure described in the following subsection. 

3.3 Integration with Respect to 9 

Consider Eq. (3.16) as mentioned above, N*(?) is a non- 
negative polynomial of fourth degree. Hence, it can always 
be decomposed as 

where are the roots of the polynomial N^, and n^ is the 
coefficient of , namely j £ * . Then the rational 

function M^/N^ can be separated into the sum of four terms 

M.ffl _ _C, c, Cc, 

~ 7 - 7 , * 7 - 7 , * 7 - 7 * 

(3.21) 

where C. are constants. Hence, the evaluation of c,, is 
1 hk 

(* ) . 

It may be noted that the problem could be avoided by 

integrating with respect to 


first. 
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reduced to the evaluation of integrals of the type 



(3.22) 

with 

V; - ^ vf <■ 2 (3 ?, + f 

3.4 Particular Cases 

A few special cases need a special treatment or can be 
evaluated in a simpler manner. They are considered in the 
following Subsections. 

3.4.1 = 0 in the integration interval 

As mentioned above, is a nonnegative polinoraial . Hence, 
if NM/O = O with r ) t real and within the interval (-1, 1), 
then has at least a double root at . In this case, the 

decomposition in fractions is different and Eq. (3.21) must be 
replaced by 

_ c, __ c. c, q 

'<■<)- rjf * 7 - 9s * n-% (3 ' 23) 

However, as mentioned above, when this occurs, \3. e ^3.,J~0 and 
the numerator is also equal to zero (single root) . This implies 
that C£ in Eq. C3.23) is equal to zero. Hence, only integrals 
of the type given in Eq. (3.22) are involved. 
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3.4.2 N 4 h 0 

The case = 0 is also possible. This implies that 3-, 
is always parallel to 3, . This, on the other hand, implies 

— /j- \ 

that all the lines "5 = const converge into the point p 
This means that the element is a planar element with two edges 
passing through the point p 1 ^ , which lies outside the element. 
In this case, the coefficient C^. is identically equal to zero, 
as evident from Eq. (1.16). 

3.4.3 Trapezoidal planar element 

If two edges of the element are parallel (trapezoidal 
planar element) the integration can be performed in a simple 
fashion. Choose ^ and such that the parallel edges are 
the edges 0 - t i or 

%■ ■- (f.i (P,t ?,) W.24) 

In order for the edge to be parallel, the vector p ( and p 3 
must be parallel, namely 

p, = y a 

h • 


(3.25) 


(where y = 


?, 


If, 


is a unit vector) and 


a, --a., = (p, „ >p s ) yy) a (3 - 26) 

Substituting these into Eq. (3.15) gives (note that ) 


A.- 


% • u p. • u « p. 

I 1 1 la. * o, I 1 


(3.27) 
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integrating with respect to f) , one obtains 


1 


P 



P. • U » P, 

2 TT 


ill I ' 


(3.28) 


Noting that (.see Eq. B.l) 

HI*’ -- | | * W | + If * W | 


(3.29) 


one obtains 


i f 


2T1 |f | 


T 


(3.30) 


wi th (Eq . A . 1 ) 


T . If I l a .-t 9 ' , = < J n 


(3.31) 


If I 


-I 


/ JeTTJ r 

where it has been set 


tan 


(ea-fed)9+M-/>c) 


J n J 


wi th 


£ . u - ( p 0 + ^F,)- w *■ 7ffc + *P»)*u 
r a + i> ? 


ll.uf. i|i 2 - u-“) 

= I Po^^F<r + 2 <p.**p.WWfr>7 4 1 K+fP* 1*7 

= C 


(3.32) 


(X 


C. ‘ U 


C* » u 
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oL . 

C* . c„ 



C„ ' c* 


V = 

C> • C 4 


c * 

JL - ft* 

-[C 0 xGlj 5 l p 0 A U | 

a - 

ft - ab 

- ( C , X u) • ( Ci * (/ ) - ( p 9 x u)' ( p a X U j 

e - 

Y - ^ 

= I C x Ct S = 1 p k * u | 

f . 

l - “ 11 

5 (3.33) 

r *■ 

c. 

• p. + ? ?, 



where 


c. -- ?>■*■? ?, (3 - 34) 

The results obtained above can be rewritten in a more compact form 
by noting that (see Eq. B.l) 

,2 


e.c - d 


P* * U I l X U I [( P„ « U ) * ( P 2 A u j] 


= 1 ( p.« Ci)* (p.x^l - 1 p 6 M X 

-(p/* p.'-u) = (P»* P.-uP * f ‘ 


i . 


where the superscript N indicates the part normal to u; 

v * ?, - ?:« (p;-v a) 


•a. 'a 

Also, it may be noted that 


(3,35) 


(3.36) 


1 f 

(ad-bc)] = [&(ey+J)- b (<lyt c) J 


I 
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=. (a {X *) + ?)- b (?? + J) J yjyj 

= [(c..ii)(^?cj)'f.- Cc' - m)U. '■ycj-c.) - 

= [(c.-CiH f- d)-(d-M><f- f -)] ^|f| 

(c~ <cj-(f»u) _ _ (f 

ir fc l f • w« p,| r h l J • w * c»l 
_ _ (f«cj-(w.) _ _ Oj^Ji (MJ 

rjf-i,«c,l JFT If •«.*«.! 

(3.37) 

Combining Eqs. (3.30), (3.31), (3.35) and (3.37) yields 

Id = — ~ sitjn ( | . Q, K Q 4 ) J (3.38) 

with 

J = + torn ' 1 - ( f xa ~Hf»^ 

1 I I I % • «. * M 


(3.39) 
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3.5 Hyperboloidal Element 

In this section, it will be shown by differentiation that 
the result obtained above is valid for any hyperboloidal element. 
Note first that 


% 

*9 

d 

Ty 


- 4* 

- P, 

« o 


Next consider ( + 1 indicates the sign of £ • fl, x Q* ) 


(3.40) 


s J j> -fcbY ) 1 — ~ ^ * lal 

n " *9 m (- % •«.**.) 


; i 


I + filyillliliiLl 1 


'Tf-J (j. Mty 


x. | [(«,< «,)• Cf »«.)•<• (|«. (j« «,)•(■«.* «.)] — 


111 ( ^ • a, « a.) 


, (fri, )■(}**.) 


~t- jr-4M>+ Jf7| («,. 5, X «. 1- 1 • ft , tfl 
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_ - r* (j - a, « I 

“ + y‘( }.«,*«,)* - »*,)* 

x |((a J < a ,)•($.* pj-f £*«!)]( £•?■)(£ •<?,*«>) 

- ( £•£««,)♦£•£ ( £■ ?>‘«>))} 

= ; i | 

*+ ( t&J) 2, r 

x | ((. », < «,)• ( %* a J (£• D - <p «'X y «•>< «* ■ • £)) 1 1 ■ «• * * J 

*Ut*hni*wi--«’ , *3- ( P *•* P ^ P P i 


( 3 . 41 ) 
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Next note that (see Eq. B.l) 


fj- <p.a l *aS*C(i*V-(i*aJ]‘ s I T I l-M-l* < 3 - 42 ' 

r 

«, r I £ * a. r - ( (£• P ( £ • «' « (( V ijj 


= ( (f 4>( V «, ) - ( i- a*/ ] ( ( f • P l «* • «0 - f f • «0 3 

= ( f. f/(a,- $,)(«>■ a,)^{\- QpJ4^S- 


_ ( %■ p) (f . a, * sj- (£• £ M • ?.)- ipaS^-) * xt-Vti ’ i'XP-iXP 

* (£•£) J - (£•»■*«*) 


- (GrQjipG.Xf-^ 

— ( £. j.) { ( £ £) ( «, * «»> rt' *«)) - f - • f ■ «- x «•) £ • ( « i * %>) 

- e, (%•(>,)- a J%- a,) | } 

= (£• !•)[ | W4.f/- 


( 3 . 43 ) 
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since 

A * (&>* c) = fr ( A ■ c) - C (A ' fr) ( 3 . 44 ) 

It is worth noting that, as shown above, 

(pQ,*v * 

( 3 . 45 ) 

Next note that (see Eq. B.l) 

?.)• i U«J(p V- ( f* «■)*■ p ( f ■ ■ «, * <?o 

+ «.)-($-* p h* L)j (f-f) 

— («>• $0d -p)(pf) 

- 1)( P «,)} ( pj {p«.x «>) 

+ {Uv ~«j-(p p(f f,j\ (%■*■' v 

- ((p %) <«,* ej-(p ij(p s.JJ (i- f,x gj (|\ 

= !•)(«,■ (£• «0j ( p a,)( y a, * &) 

+ {(p %)((*>• b)<-p $.*&)-(«,' i J( pK*K) ) 

- (£• %)((%- p £• «,* «V- (f s,x pf» l))k. I 
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= -1^‘a.r (f.- «,)(£•«.*<?.) 

(3 _ 46 , 

since (note the change of the order in the triple product) 

(%■' %){&■?>)( «,)-(*,■«>) ($:*«.. p 5/ )J 

- ( £ • <?J (( \ • p» ) ( «- ■ «)- < V «<K% ’ «" • f*)J 

- ( r «*) Of- «,;J 

= £’ P)««, 

= If* «X%- ?>*«' 

(3.47) 

Finally, by combining Eqs . (3.38), (3.39), (3.41), (3.42) 

and (3.46) one obtains 



28 


— — i<kr\ 

i// 


/TT (*-i- «.*«*) 


i %>. «, r ! i * «, r &' a ' ,( v a ~ ' * 5j) - ■ i * p, >J 


r i f ^ r 


(3.48) 


and, according to Eq, (3.38) and (3.15) 

-isj -i i 


2 Ep 


2IT 3 


-2-rr A i> 


2TT K 


if *«.i‘ 


(3.49) 

in agreement with Eq. (3.28). This completes the proof that 
Eq. (3.38) is valid for any hyperboloidal element. However, 
for the sake of completeness, the derivative with respect to 
is also performed (note that -2-~- = q, . ? , ^ * — - p and 




0) 

2 TT 




® In 






Jyt 


~ *5 V jyr i ?. o i 1 > 


ifx«,n r#:. 


' ) H I M.( 

^ f ' 3 •{« &)- r }■ «.)( f • «. * «3)J 
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+ 


2 (tf 





ffTjr 1 -' /T* r<^‘h) ■ r ^-p( a ' ‘ ^ K H) 

-(a, • a, )( f • a, « «J- (£• *,)< «, • «, x «0 -(§■ •* 








£ * *.* 


in agreement with Eq. (1.16) 


( 3 . 50 ) 
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SECTION 4 
SOURCE INTEGRAL 


1. Introduction 

Consider the evaluation of b hk - According to Eg. (2.21), 


b, , can be evaluated as 
hk 




• StI ( d ? 


(4.1) 


where 


aij_b 


iH + c (p) 


(4.2) 


with 


a = - i.iJ - p, x p, . Z x 

-- (.j, *. -b, i») + ? - (p, *ph)* h'A 

X - - 1 

f = 2. -i, = (f >+9f>)-(f.~9K) 

Y • a. •£. = (p^9h)-(?.*9h) 


(4.3) 


Note that Eq. (4.2) yields (Eqs . A. 9 and A. 10) 

A A . . 


A 


A / 

a 


5 f 




* / fT? 








(4.4) 
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The integral obtained combining Eqs . (4.1) and (4.4) can be 

easily evaluated numerically. Note that b is a linear function 
of YJ while ^ and ^ quadratic function of ^ . However, 

for elements with two parallel edges, the integral can be 
evaluated analytically, as shown in the next Subsection. 


4.2 Trapezoidal Planar Element 

Consider an element with two parallel edges (trapezoidal 
planar element) . Choose ^ and such that the parallel 
edges correspond to ^ = 11. For this case 

p, = 4 

p, * y w <4 ' 5) 

- V 

(where U s — 4r-r ) and 

Ip, I 

2, * Cp, ^n>) • (4.6) 

These imply 

A 

a - o 

X - (x + wf 

£ * {X + • £.) 

t * k'J. (4.7) 


and Eq. (4.4) reduces to 

A 

As * j== ^ ( JT V" + if ♦ f) + C (y) 
+ C(J) 


(4.8) 
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Note that 


Ip* I < Ip, I 


(4.9) 


otherwise the boundaries = + 1 would be crossing. This 


implies 


^ > ° 


(4.10) 


and hence, Eq. (4.8) reduces to 


- JU [ f t + + M - t C,(j) 

= & + (<!» + ■/ il, ) * u ] + Ci (y) 

= & M I + ccj; 


(4.11) 


C,( 7 ) = c ^ 


(4.12) 


Consider 


Sf^A-sdi? » ^ [i" (I I- 1 + V') ir ) * c ll >) 


(4.13) 


ciD-bfC'Wi f c 


(4.14) 


Integrating by parts yields 






« (?-*)*. < t|! ♦ a)- 1 * c r„ 


(4.15) 
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where is a suitable constant. The choice for ^#-is 

discussed later on (see Eq. 4-. 30). Note that it is possible 
to write 

( 9 - 9 t )(K+*r,) - ( 4 .i 6 ) 

wi th 

h = F. * ? « + f % * sf, ) (4.i7) 

By using Eq. (4.16), Eq, (4.15) reduces to 


■si 


5 ~ 9+) ( { %-j + p • ft) 

_ jufij gjj) I , f- 
/ Ifhf'M ifi ' / f * ifl + f-a 




+c (?) 


tiJRA 1 


with 


ffl 




(4.18) 




(4.19) 


Next note that 


S1 ) ** l|. 1 7 j-u | Ifl ^ ^ 


1 if i-f.* 


) 




%*■% - (fr,'U)(£ j i) j n , -(f,- !)(£•«) 

lfr-(f-ar > j Ifr-(j .«)* Ijr c >^‘ 4 - 2t » 
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Applying the general formula (see Appendix B) 

C x p) * (#*C )(&'$)- c) (4.21) 

one obtains (note that VL • U = 1) 

lf-1 - I j.. M I * | f * U I - |£. ' ^ I 
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since the integrand in Eq. (4.26) is independent of ^ . 
Finally, note that it is possible to write 


h - - f," * %l m 

(4.27) 

where is normal to both {/ and P 3 or 

&/• U ~ ° , + h) (4.28) 

Combining Eqs. (4.17), (4.27), and (4.28) yields 

(<p. vp.*-? f>)~ fjui ■ w »» 


(<P.*?pJ - 9* ll hJ- (p>+fh)*° 

(4.29) 

By using Eq. (3.33) and noting that w • Of - l y Eq, (4.29) may 
be rewritten as „ „ , 

&■ + %t- u ^ 0 


f - %f-K b *o 


that is, a system of two equations for the two unknowns 



which has the solution 


?*- and 



f - ab 

r- b* 


£ 

e. 



aX - bp a(-bp 

K - b* e. 


(4.30) 
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Finally, combining Eqs. (4.25) and (4.27) yields 


isi - — rrr?:T— :r, r h ttt 7 + 


I % • M I 


Ifl 


*hj Iff c * ( ? ) 


(4.31) 


= -r- p. 


f ? • m a 9 

. -a- < 

? T j 

r j 

)/£*«/* l|-l 

t 1 


1 1 



since . y -q (Eq. 4.28) 


&.*• £ * h"- p. 


(4.33) 


For is normal to u and p t + p^ (Eq. 4.28) . Note that 


(see Eq. A. 10 ) 


^ - j ( dy 

= Jf ^ I ?■ I +■ 


(4.34) 


with ^ and jjj'' given by Eq. (3.33), while the first 

integral in Eq. (4.33) is proportional to Ip given by Eq. (3.28) 
Finally combining Eqs. (4.18), (4.32), (4.34), one obtains 

Ijtf-W = + -j==4(/f \\\*Xn*s) 


~ ?'• P»777 JJ 


(4.35) 
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with J given by Eq. (3.39). The results obtained above may 
be rewritten in a more expressive form. Note that 


7 - 7, 


4 £- 7 t d ( c>xti-c i xCi}y + (c 0 xii)- (c>.xCi) 

^ & (C v xm)*(^>* m) 


» (4 X X4,) 

iCj.* a j x 


(4.36) 


af - bp _ ( c >‘ u)(Cj r cj- (Csu)(c,-cJ 

*£- I Cl A & I* 

_ &*C>)'(Ci*Ci) _ 3 ,) I 4, 1 

I X W ) I K Q x | 


(4.37) 


and finally 


fv • P. 
Ifl 


?/• r 


I Pc • n | 


i f i l m * pj 


1 ?*' n |* 


I P» * n 1 

I P**n 1 1 m * F*| I w a p 4 | 


[gif 1 %• • g, x 

15, x 5a r 


(4.38) 

Combining Eqs. (4.35) through (4.38), adding for convenience 
the function of 1 - , and noting that (Eq. 4.7) 

b ~ M * £ * /L* (4.39) 
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yields 


I, - & 


2-TT 


^ (/t + yj * 5, ] *+ 

jKmJT ■'.{.fj-ttij} 


5 > 1 ^ (.(JTiJTi- 

+ ( £*«*)• (v^) Jf*! * 

H ^ j. * 

- 1 «,* «J j ) 


(4.40) 


4.3 Quadrilateral Planar Element 

In this section, it will be shown, by differentiation, 
that the results obtained above are valid for any quadrilateral 
planar element. For this element, the normal 

- 5: , * 

n = -j-z — zz— 

i Qt * I 

is independent of ^ and ^ and Eq. (4.40) reduces to 


(4.41) 
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I s = _L n-X *«,-« ~j <U\ (|«,||f| +|-«,) 

■+ £ x -j4-j 2* ( |7, 1 1 1| + |3«») 

- If • " I J j 


Note that 


2> 

*5 


* O 


and 


r~r ( fc* Mixtfrn * * 


M |5ill jM *£•«,) 

z> -af v 


1 f *■•*' /? I , .0 .? 1 - Jv a ' 

ir,iifi+i-*, ijfj 1 '' ' 'J 


Hence 


•Z> 


f «/ 


( v?> 77] ^ i|-) + $-* 


w - "' 4< ’ 


(4.42) 


(4.43) 


(4.44) 


r* 


(4.45) 
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Similarly 


f C \ c? A • n — — bn ( I Ai. I j \ \ i- 3.' 5 a ^j 

L r |«,| r ^ 7> 


„ ^ / jrj l£uZL) = ^ 1 n ^ f*F>-n j - - 

Jfj f Jfi fR 




r 


Furthermore, noting that 


(4.46) 




"3 


n a-o 


and using Eg. (3.50) yields 




■ilf-nU) • 


(4. 47) 


- M 


(4.48) 


Finally, combining Eqs . (4.42), (4.45), (4.46) and (4.48) 

yields 

** ^ . iii^L t <pi*h)L -A. + 

na, ( r 3 

^ f • 4 f (bn) %'4\xi x 

+ ^ if.'*,*") ^ -pi ; y; 
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(4.49) 

in agreement with Eq. (1.15). 

4 . 4 Comments 

It may be noted that triangular elements are the limiting 
case of quadrilateral planar elements. Hence, the above for- 
mulation is exact for triangular elements. On the other hand, 
if the derivatives of n with respect to ^ and ^ are negli- 
gible, then Eq. (4.42) can still be used. It may be noted 
that the error introduced by assuming n s o within 

the element is of the same order of magnitude as the one in- 
troduced by using constant-potential elements. Smaller ele- 
ments are thus required where 0 is varying rapidly, that is 
at the leading edge and tip, where, incidentally is also 
varying rapidly. 
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SECTION 5 
THE WAKE 


5.1 Dynamics of the Wake 

As mentioned in Section 1, the surface in Eq. 1.1, 
surrounds the body and the wake. The effect of the wake, 
disregarded in Section 1, yields an additional term in Eq. 
(1.4) , given by ^ ’ ^ 

-p-da- (5-D 

<r 1 1 

with 

( 5 . 2 ) 

This represents a distribution of doublets with intensity 

. The geometry of the wake is not known. An iterative 
procedure can be used to solve the problem: consider the 

surface of the wake divided into small elements. Assume 
initially that the wake is composed of straight vortex lines 
(see next subsection) ; then find the values of and then 

evaluate the velocity at the corner of the elements. Find 
a new location for the corner of the element such that the 
elements approximate the stream-surface emanating from the 
trailing edge and repeat the procedure mentioned above. 
Needless to say, convergence of the iterative procedure should 
be verified, numerically, if not theoretically. 

Finally, the values of a ^ at the centroids of the ele- 
ments can be obtained as follows. Consider the Bernoulli 
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theorem for potential flow 




dt 2 


Since no pressure difference can exist across the wake, then 




($«-$,) ^ ^((d <iv *$.)-( *4y *<!>,)] = 0 (s-5) 


This can be rewritten as 


-*• t (5 (^+ <h)* v $,)] - o 


(5.6) 


A <§ “ O 


(5.7) 


where 


— t Q * ^3 
3-b •" v 


(5.8) 


is the total time derivative obtained by following a particle 
having the mean velocity, , 


v <|> w + v <k 


(5.9) 


For steady state flow, Eq. (5.7) reduces to 


^ <M = o 



44 - 


7 

(in agreement with the result obtained by Mangier and Smith ) 
or 

^ - 4 * * ***** ( 5 . 11 ) 

along a streamline. Hence it is convenient to use elements 
with edges approximately coincident with stream lines. Then 
the value of is the same along the strip obtained by con- 

tinuing all the elements between two streamlines. This yields 
a simpler set of equations since only the value of <&§ at the 
trailing edge (rather than the values of at the centroid 

of each element) would be involved. 

It may be noted that the above derivation is exact in 
the sense that no small perturbation hypothesis has been made. 

It may also be interesting to interpret these results in terms 
of velocity: the vortices are parallel to the streamlines, 

the total vorticity between two streamlines (equal to the 
difference in a <j> between the streamlines) is constant, while 
the intensity of the vortices decreases if the vortex lines 
diverge. Note that the vorticity is given by 

^ 4 \) (5.i2) 

where s is the arclength in the direction normal to the lines 
- const. 

It may be worth noting that (Ref. 1) 

I! ^r(-r) ds = n (5 - 13) 

<sr 

where XI is the solid angle of the surface as seen from the 
control point. Hence, the important factor is the contour 
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of the strip (rather than the shape of the strip) since the 
solid angle depends only upon the contour. Note that the 
doublet integral is exact for any hyperboloidal element. 

Furthermore, it should be noted that after a few span- 
lengths, the wake-sheet rolls up into two vortices. Hence, 
after a few spanlengths , the strips can be replaced by two 
concentrated vortices. 

Finally, it should be mentioned that, as shown by 
7 

Mangier and Smith ,"the vortex sheet shed from the trailing 
edge of a lifting wing with non-zero angle, in inviscid 
subsonic flow, leaves the trailing edge tangentially to the 
upper or lower surface, or exceptionally, in an intermediate 
direction. The exceptional intermediate direction is possible 
in three circumstances only: either there is no mean flow or 

no shed verticity at the point" or both. If the shed verticity 
is positive (negative) and the mean flow is outboard (inboard) , 
the sheet is tangential to the upper surface, otherwise is 
tangential to the lower surface. Note that is has been impli- 
citly assumed that the flow leaves the aircraft at a sharp 
trailing edge. Shedding of vortices from the body requires 
the use of viscous flow equations (Ref. 1) and is not con- 
sidered here. 

5.2 Simplified Treatment of the Wake 

A simplified treatment of the wake (used in Refs. 1-5) 
consists of assuming that the wake is composed of straight 
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vortex-lines emanating from the trailing edge and parallel to 
the x-axis (direction of the flow). For this case, the 
surface of the wake is divided into infinitely long elements, 

A 

<S1 , with two edges parallel to the x-axis. These elements 
are the continuation of the elements of the wing having an 
edge in contact with the trailing edge (Fig. 4) . 

Hence, by assuming that the value of a ^. e can be approxi- 

A 

mated by the value at the centroid of the element S' the con- 
tribution I w (see Eq. 5.1) is given by 




(5.14) 


wi th 


• i* - i iWf * ("b) ( 

for the elements with an edge in contact with the trailing 
edge , and 


= O 


(5.16) 


for the others. 

In order to evaluate the integral in Eq. (5.15), it is 

A 

convenient to consider that the element <5^ is the limit of 
the parallelepipedal element obtained by truncating the 

A 

element <5^ at the finite distance (Fig. 5). The limit is 
obtained by letting 

f , • / * * a x (5.17) 


go to infinity; note that u = i since two edges are parallel 

to the x-axis. Note that (see Fig. 5) 

O _ V P+ P- _ 'p 

lo 11“ 2 “lm 

P, * f 1 -5, 
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P* - P- 


2 


0 




(5.19) 


It is convenient to separate the contribution from the 


trailing edge ( = - ) ) and the edge that goes to infinity 

( = i ) : 




where 

(note that 

£ = ? n ^»i)t+ r ) h 

) 


S a 

( |. • «, x §J - 

' a y rv <£/*< ft) 

with 

P«j = 

p* - 7 Pi 


while 

(note that 


) 


J w (. t,<)) * 

x-»- Mfi 


, A* 

\ I ?- 1 (£-2 

^ L v -/ - i&+nd^l-tip„i-2xi),fj,] 

. r 1 - (Pm A !)■(! » Pj) 

I hi 

fed - k 

I f„j ■ <- * P J 


(5.23) 
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and similarly (note that (-£ ^ = _ j? 


md 


JwK?; 


Ifl J 


■5 


^UA ' 


r "fotH 


(£*X) » (f * Rj) 

l If I I f'£*Pdl J f 

-chid**)' (fi* * ft) i 

I find I J fmi ’d'*?A >■% 


= ~t&r\ 


~'(P*td'Pm<t')( A '‘ft)~ (pmj' Pd )( Pt»s * **- ) 1 
I Pmd I l Pd | 


(5.24) 
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SECTION 6 

COMMENTS AND SUMMARY 

6 . 1 Comments 

As shown in the preceding sections , the use of hyperbo- 
loidal elements is quite cumbersome unless two edges of the 
element have two parallel edges. In this case, the element 
reduces to the trapezoidal planar element considered in Ref. 

1 tin particular triangular) . 

Hence, it is convenient to divide the surface into tri- 
angular elements {which are a particular case of trapezoidal 
elements. For flat surfaces {the surface of a building for 
instance) quadrilateral elements are more convenient. 

It should be noted that the doublet integral is equal 
to the solid angle (multiplied by ) of the surface <5^ 

as seen by the control point . Hence, the correct shape 

of the solid angle depends only upon the contour of the sur- 
face <5 |t but not upon its actual shape. This implies also 
that it is important to use triangular elements, rather than 
tangent plane approximations to the hyperboloidal elements, 
which would yield discontinuities in the surface and hence, 
the total solid angle would be changed. A measure of the 
accuracy of the method in evaluating the coefficients 
(Ref. 1) is that the sum of the coefficients should be equal 

A 

to -1 (solid angle multiplied by ) . Note that if the two 

2TT 

hyperboloidal elements in Fig. 3a are replaced by the two 
hyperboloidal elements (in the limit triangular elements) in 
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Fig. 3b, the centroids of the elements move very slightly. 

This is convenient for the evaluation of the pressure coeffi- 
cient Co - _ 2 by central difference (since centroids 

rx 

of the elements lie on the same wing section if the elements 
are bounded by wing sections). On the other hand, hyperbo- 
loidal elements can be used if the derivatives of n with respect 
to and are negligible (Section 4.4) in the source inte- 
gral (the doublet integral is exact for any hyperboloidal 
element) . 


6 . 2 Summary 

Assume that the surface has been divided into quadrilateral 
elements. The elements are described by the vectors p e ^ p , 
from which one obtains 

f F* + fi+ 9 + hf- F* U) 

a ' = ft + 7 ft)** , + 

Then the approximate solution of Eq. (1.1) is obtained by 
solving the linear system of equations 


with 

^ Lk ~ Qk " 
- I k K 


( 6 . 1 ) 


( 6 . 2 ) 

(6,3) 


where w ik represents the contribution of the wake discussed 
in Section 5, while and b^ k were derived in Sections 3 and 4, 
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For convenience, 

coefficients C., 
ik 


the results are summarized here, 
and are given by 


The 


biK - 1. (1,1)- 1, (1,-1)- XsC-1,1) . I. (-1,-1) (6-5) 

In Eq. (6.4), I is given by 

Ip = - r~ pfn ( y . Q,a «J J (6.6) 


with 

T - Q>) ‘ x $i.) 

Ifl 


(6.7) 


On the other hand, in Eq. (6.5), I g is given by 

. Is ^,l) = -n-l l s l^,n) 


( 6 . 8 ) 


with 

j - C • n) j== In (Jit, ■*,)(£■$-) t pS, 

+ ^ ,Jf?| ^(/J + 

- I F’"l J) 

(6.9) 

Equation (6.9) is exact for quadrilateral planar elements 

and may be used for hyperboloidal elements if the derivatives 

z> n p n , , 

, — are negligible ; this is possible if small 

' Tt rj 
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elements are used where the surface curvature is high (leading 
edge and tip) . 

Finally, the coefficients w^ in Eq. (6.2) are given by 

Wi K - D (6.10) 

for the elements without any edge on the trailing edge. For 
the elements in contact with the trailing edge, assuming that 
the wake is composed for straight vortex lines, the coeffi- 
cients w^ are given by 

with 


a 9 } 

I Fm4 '** frl 


( 6 . 12 ) 


and 


Jw(-t - 9 ) = 


(6.13) 


where i is the unit vector in the flow direction, 

5 m P. 1 P- 

\ Tf\ O 


P * * ^2 
? '-i * P* 7 h 


(6.14) 

(6.15) 

(6.16) 
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where p + and p_ are the location of the trailing-edge corner 
points p r ^ of the elements with respect to the control 

point p fK) 

p. ‘ P™ - p" 1 ’ , p. * p TS "_ p w (6 - 17) 


6 . 3 Exterior Neuman Problem 

The results obtained above can be applied for the solution 
of the Laplace equation outside the surface <S , with pres- 
cribed normal derivative on ^ and usual regularity conditions 
at infinity. In this case, Eq. (1.6) is replaced by 


* - 


* 




<1 <$. 


d- I [( n>v(y)6<s K h*i,2. 


and correspondingly, Eq. (6.1) is replaced by 




Y\ 1 
k ) 


(6.18) 


(6.19) 


with 

wake) 


a ik given b Y E( 3 • 
and 


( 6 . 2 ) 


(with 


0 if there is no 


b.v = I.d.l)- I, (1,-lJ- i s M.1) + (&. 20 ) 

with I s H, *?) given by Eq. (6.9). Clearly, if the boundary 
condition is given by Eq . (1.3), Eq . (6.19) reduces to Eq. (6.1). 
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APPENDIX A 

USEFUL INDEFINITE INTEGRALS 


A. 1 Introduction 

For the sake of completeness, some indefinite integrals 
used in this report are verified in this Appendix. 


A. 2 Doublet Integral 

For the evaluation of the doublet integral, I D , the use 
of the following integral is required 


with 


J - 


cl -t-ky 
<-■* z.dy+z.y*' f 


-d 


7 


tv I 


l 

4 


'C *■ £ (a-*- by) +■ ( c+- i'dy-b e ;*;h 
= t («*c ) r 4)p*( &)p* ] 

Differentiating Eq. (A.l), one obtains 


£ 


(A.l) 


(A. 2) 


f _d_ -fa' * ( (ae-btt)9 i ■((t.i'bc) I | 

Jec- d*’ 6*} o J&c- F ? ' 

j 1 {(&e-U)?+ M-bc)} j(b+i) r )+ (*b+ 6)) 

/itc-d* i + J ec-r 

l rite- S' f > 

(e.c-A)X * £(ac ~i>d)9 + (a.4 - l>c)J * ( 

A -r l>9 I 

c ■ +* ^ 


(A. 3) 
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since (ae.- 

s (ae-bij {(&b + 4)y +■(* + *■)] - (*.&-U)[0?+t)y t (*b-ni)] 

® 4)- (*4- be) (b+e)} ^ + ((«e-- b4J(atc)- (*d-bc)(*b+d-)] 

s j a. de - a b*d - M *- oC’^d -a4jz + b* c- + bee J p 
+ t a ce - - i>c d - a bd - « <4 + ab r be 6 J 

- Cby + G)(Qe-~d*'i- b*c- + c c - *-a.bd ) (A. 4) 

and 

(ce^ — dpY + ((&e~bJ)j?‘h(Ael~bc)J 

* ((cz-r)U>*t*)+(«e-M/ ]p\ ((ce -A*)(&b*d) + 

(ae. - bd) C&d - be) ) +. C (cz- dVi&’+c) + (Ad-bt ) J 

= [ b*c.e rc e*- b t d*'- 4 « + Q^e.* -iA.bd*. + b^d* J J 

+ £ Q-bce.* c d& - & bd — d ^ -4-^ d£~~ a i?c e. - G b d + b & d ] *■*) 

t [a*ce<- t J e - *-ebc d + b c J 

= (e}7*+ *-dpr&J(& C- d*+ be +• Ce-zabd.) Ca,5) 


A, 3 Recurrent Formula 


The recurrent formula 

^t T *9) n 'yfy •({+*?) r_ 


n - 1 




*-»• 


y 7 


o 


(A. 6) 
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with * s J A +*?) + if 

is easily verified. by differentiation 

.n-l 


(A. 7) 


h, (?*»?;“ 7 TTZ^TTf 

- V.n- ff)" J J- f?‘ <■ (?'■(?)'"' .- I ** l 7 

(p + nr* , 

‘ * oV * «»♦ < 7 / J 

s J*»v!rr (* 1 w** j 

- »i ~r r (*-1 K&X-^fp + tp)* *~L 


In particular, for n = 1, one obtains 

( 


[Itl+ty+ry ' * ‘ 7 


(A. 8) 


(A. 9) 


Note that for n = 0, Eq. (/A. 6) is not valid and is replaced 
by 


’ | | 

T"*? W*'* r ? + <^ 


(A. 10) 


For 


j=4^ r ^^ 


/F J7 ^ r?*p 




r 


(A. 11) 
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\ 

APPENDIX B 

A VECTOR CALCULUS FORMULA 

A convenient formula used extensively in this report is 

(a- ih'c)[p>§) * {ft* f?)'(e>x i) (b.d 

which is easily verified in the following: 

{A • S)(c- t>) - (ft* c)(b-v) 

= (ft*e>K ? fty^ + A $ 3s)(C*P K + C } Vf 
~ * AyCy + A % C f )( e > * P* + ^ P*) 

- p x - /fy *»• + Ax a* 

+ A* 0x Cj P$ + Ay C x Px i- Ay pj C„ p x + Ay &y Cj Pj 

r P* h V }~ ^ A*pi £' Px ♦ ^ 

t Ax C* Py P^ + Ax <r x pj ^ Ay t3»x P* 4 - A} Cj B>x Px 

” & > P >- ^ B,P # j 

(&• 

(a * p ) ♦ ( I X c ) 

fl j * I / * I * \ 

Ax Ay aJ. #x J » 

‘ P« Pjf P>/ \ ^x / 


while 
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( A) - A) Vy ) C C * - ) 

+ ( fyv* - />* v } )(b } c*~ &*c s ) 

* (fixVf- fry Vk) (&* " &f c *) 

Ai - *} b , c , p ,- (*s 

+ P* + Ax &x C ^ Pj — (/fy 

+ ft* fix Cy Vy -t- Cj, P* _ (fyy 


e>}C}l>i-A$ 3y Cy Of) 

c ^ 3 x P* ^ ^ ^ j ) 

Cy (3x 'Px +• Ci TPi^ 
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APPENDIX C 


SUBSONIC OSCILLATORY FLOW 


C.l Integral Equation 

In this Appendix, it is shown how the results obtained 
in the main body of this report can be extended to subsonic 
oscillatory flow. Introducing the variables 

X = Y- _1 Z.- JL T-.f. ft* mJL lc J 

pi l l r l />«» '• 

A 

and the complex potential (fi such that 

r ^ , ;n(r+nx)j 

§ (•.}*)* t (YXZ)e / (c .s 


the integral equation for the subsonic oscillatory flow is 
given by 


h-i 


r- ^ j' -l&f? s. - -,'Y1R f 

ii. S— + ef> JL e )]j£ 

dN K ' 9N \ R / 


where Z" surrounds body and wake. 


C.2 Boundary Condition 


The boundary condition is given by 


v S-V <P = - M. - x. 

**? / dx 


(C . 4 ) 
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or 


(7 $. V <i ?l. 4 ± &L + !¥l te-ex.o 

xrT. xVi / >T a a v/ R 2 - ax' 


A 1 dT /3 dX p 

where cj) and c j) are such that 

<£> z (7 x (j) = £/ l (X+ <f ) 


(C . 5 ) 


(C. 6) 


Next, assume that the motion of the surface consists of small 
harmonic oscillations around a rest conf iguraiton, that is 

$ = S(X.YJ) +3(X,r,Z)e' nT (C.7) 

Then, setting 

<£> - cp (X, Y, Z) + <p (x, Y / Z) e (c . 8 ) 

one obtains 



(C . 9) 

Assuming 

5„ - O(i) cc.io) 

3S„ c Ole) 

?X 


(C.ll) 
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with 


and 


KS . 

= of/J. 

(C. 12) 

w ° 


r*J 

s - 

0(£ 2 ) 

(C . 13 ) 

.a = 

0(lJ 

(C . 14 ) 

1! 

ivm*< 

0(1 J 

(C . 15) 

e that 



<£ = 

ofej 

(C. 16 ) 



(C . 17 ) 


1 2 Jl i 4 

Neglecting the terms which contain € (of order £ ) and 

separating the steady from the oscillatory terms, one obtains 


, _i ds a + m* gg° d A ^ 

Xr° ' * /3 3^ 


o 


(C . 18 ) 


(C . 19) 


-f- M z f d&o d(p ds 
ys 2 ( ^ ^ 


- o 


A 


Introducing <b such that 

^ -a iJlMX 
(ft = f) e. 


( C.20 } 


Equation (C.19) reduces to 


y /JZA7X 

VS ‘ V 4> € ^ 

WJ 0 


^ ^ ifLWX ~ „ , 

,'JLrf & e + VS -7$ 

~Z T f 7 xy?/ * 








i 

-f 


M ii ] , o 

2X ?x J 


JlMX 


(C . 21) 
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Finally, neglecting terms of order £ in Eq. (C.18) and terms 


of order £ in Eq. (C.19) , one obtains 


VS .17 i ■- 
V & .\7 4> 

m o xyg T 


1 

dX 


■ a a < ^ 

* f p -Jx 


In particular for 


£ = ±i U- *.(*•?) -*(*'$)* 


ico6 


l 


(C . 22 ) 


(C. 23) 


(where the upper [lower] sign holds on the upper [lower] surface), 
one obtains 


and 


where 


4 , -*■(*■ 


aJ 


73 - + { * (*■?) 
& 1 H > ■ 1 *' 


A 


H E <&■ W - - N (ikS + 21 ) 

dN / vs./ 2 **■ 


A 


j2 = OJ l 


n 


u. 


(C . 24) 
(C . 25 ) 
(C.26) 


(C . 27 ) 


(C . 28 ) 


Equation (C.27) gives the value of d<j>/dH to be used in Eq. 
(C . 3 ) . 


C.3 Pressure Coefficient 

The pressure coefficient is given by the linearized Bernoulli 


theorem as 
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CL = _ 2. (2/ * |« 

r u zK 9 b 

- llA M ) 

' ' (ft j?r V 2*V 

For oscillatory flow, setting 

/v 1 SIT f ,Jl(T-fMX) 


<ft = p 
T " T c 


- ft e 


^ ;a/ 


one obtains 


~ 2 (fc (ft <p + J "X 

T = ' (*? r p dx 



^ _Z ' 
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^ „ iSLhK 


-2. [c'JL ft + H 1 
j$ L At srK J 

_ . /J2X/W 

2 fi 


^ ,' 1 Z A?x 
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t nA/n \ 1 ;jz/?X 


- y 

"W 


H jL ft * 

zx \ 


Z € 

/* 


2 ^ 
ft ?>X 


2 ,siX/rt 


(C. 29) 

(C. 30) 
(C. 31) 


(C. 32) 



